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Abstract 

Dynamics of solitons of the Ablowitz-Ladik model in the presence of a random potential 
is studied. In absence of the random potential it is an integrable model and the solitons are 
stable. As a result of the random potential this stability is destroyed. In some regime, for 
short times particle-like dynamics with constant mass is found. There is another regime, 
where particle-like dynamics with varying mass takes place. In particular an effective 
potential is found. It predicts correctly changes in the direction of motion of the soliton. 
This potential is a scaling function of time and strength of the potential, leading to a 
relation between the first time when the soliton changes direction and the strength of the 
random potential. 
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I. INTRODUCTION 


Solitons are one of the most remarkable manifestations of nonlinearity. They are found for 
continuous systems for the Nonlinear Schrodinger Equation (NLSE) in one dimension (ID) jl]. 
On the lattice mobile solitons are found for the model introduced by Ablowitz and Ladik (AL) 
[l2i], while for the ordinary NLSE on a lattice, a mobile soliton is only an approximate concept. 
Both lattice and in the continuum disorder tend to affect and typically destroy solitons. In the 
present work this will be studied for a ID lattice in the framework of the AL model. 

For the continuous system, early numerical work of Bronski [ 3 J indicates that, a NLSE 
soliton becomes trapped in the random media when its kinetic energy decreases sufficiently 
and is of comparable size to the background potential. Akkermans et.al. 4:] numerically show 
that a soliton bounces back and forth between high potential barriers in the attractive Bose- 
Einstein condensates in the framework of the Gross-Pitaevskii Equation with strong disorder. 
In addition, like the Anderson localization of linear waves in random media, someauthors relate 
the localization of solitons in disordered environment to Anderson localization 
In absence of a random potential the AL model is 


.dip* 


dt 


= -(VL-l + 1pn+ l)(l + IV’nH , 


( 1 ) 


where ^ n is the wavefunction on site n at time t. The integrability is manifested by the existence 
of a mobile soliton solution [7, 0 ] 

sinh(p) 


4>n(t) = 


■ exp [ik(n — x) + ia] , 


cosh [/i(n — x)] 

where the time-dependent parameters x and a can be expressed as 
sinh(/i) 


x = 2- 




■ sin(fc) , 


k 


a = 2[cosh(/i) cos(k) H— sinli(p) sin(fc)] 

l l 


( 2 ) 

(3a) 

(3b) 


From (J2J) we see that characterizes the width of the soliton and x is its center. On the other 
hand, the AL equation has two conserved quantities, the first of which can be defined as the 
mass of the soliton solution 


M s = 


OO 

E 

n =—oo 


lll(l + \^n\ 2 ) , 


(4) 


while the second can be defined as the momentum of the motion 


P = i (4C+1 - V’nV’n+l) , 


(5) 


where * denotes the complex conjugation. For the soliton solution (J2]) we can calculate that 
(also can see Appendix of 0 ) 


M s = 2 p , 


( 6 ) 
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and 


P = M s x = 4 sinh(/i) sin (k) 


(7) 


Therefore M s can be indeed considered as the mass of the soliton. 

In the present work we will study solitons for the AL model with a random potential defined 
by 


dilj r 


dt 


(VVi—1 T t/Vt+l)(l T |VVi| ) T ^n'Pn 


where e n are independent random variables uniformly distributed in the interval 


10 


W W] 

' 2 ’ 2 J 


( 8 ) 


For 


1 ] have addressed t he p ropagation of solitons in 


3 Ei 


and Gamier 


13], they show two regimes 


the continuous version, two early reviews 
disordered systems; in the works by Bronski 
for the NLSE soliton propagation. In one regime, the soliton mass decays while its velocity 
approaches a constant; in the other regime, the soliton mass approaches a constant while 
its velocity decays very slowly. Gamier further applied a perturbation theory of the inverse 
scattering transform to confirm that similar two regimes are found for the AL solitons with on¬ 
site random potential (in the limit of zero randomness) [l^] . Which regime is relevant depends 
on the value of the initial mass of the AL soliton. For large /i, the mass approaches a constant, 
while for small /i, the velocity approaches a constant. 

However, we will show numerically, for the weak randomness, the large soliton will be 
trapped before its velocity decreases to zero. Additionally, we will find a regime in which the 
AL soliton has possibility to be accelerated on the average by the randomness. We will also 
characterize the regime where soliton can be trapped by the disorder using a particle approach. 


II. SOLITON PROPAGATION IN A DISORDERED LATTICE 


In this section we will study the solution of Eq. ((S]) numerically and semianalytically in 
order to develop an intuitive picture of the soliton dynamics. The initial soliton is the one 
found for a chain without disorder, given by ([2]) with 

x = 0, x > 0 , (9a) 

M s (t = 0) = 2 fi . (9b) 

The numerical solution is obtained propagating the soliton by Eq. (J8]). To save computer 
resources we use a coordinate system moving with the center of mass of the soliton, consisting 
of N sites, centered on soliton. The computation is performed using a 4th-order Runge-Kutta 
type algorithm in time, and an absorbing-wave boundary condition on sites. We define the 
following quantities: 
soliton mass 

N 

Mf) = £hi(l + |^| 2 ) ; (10) 
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center of mass coordinate 

N 

x (N) = y^nln(l + |-0„| 2 )/M s (jv) , (11) 

n 

and the second moment 

N 

m 2 N) = ~ lll ( 1 + IV'nl 2 ) , (12) 

n 

while the soliton velocity is 

v = Ax^/At , (13) 

where Ax^ is the change of during the time interval At (here we use At = 0.001). In 
addition, two parameters, that characterize the soliton are introduced in simulation: one is the 
amplitude of the soliton 

A s = max |-0n| 2 , (14) 

n 

the other one is the soliton width, defined as the minimum N w satisfying 

(N w —l)/2 

Y ln(l + |V> n | 2 )/M a W > 1 - 5 , (15) 

-{Nw- 1)/2 

here 6 = 0.01 [note that, in simulation we first find the peak position of the soliton, then (TT5T) 
is calculated with this position as the center]. 

There are basically 3 regimes characterized by the initial value of /i: 


(41) 

|i>l ; 

(16a) 

(B) 

h ~ i; 

(16b) 

(C) 

At < 1 . 

(16c) 


A. The regime ji^> 1 


We choose an AL soliton with /i = 3, which has more than 99% mass concentrating in 3 
lattice sites (N w = 3). The reason for picking up a soliton of such narrow width is based on the 
fact that it is compact enough to admit low level of mass radiation resulting of randomness. 
Such low-level radiation is necessary for observing possible soliton acceleration in our numerical 
simulation. The initial velocity of the soliton is chosen as x(t = 0) = 1, and one realization of 
the random potential with W = 0.1 is used (note that 0.04 < W <0.1 can be seen as the weak 
randomness in our discussion, below this interval is considered as the very weak randomness 
where the soliton dynamics may approach the one in the limit of zero randomness as in |l3|). 

Assuming the random potential is a perturbation, the approximate equations for the various 


parameters in this potential can be derived following the work of Cai et.al. [7j. The resulting 
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equations derived in Appendix A are 


(17a) 

(17b) 


/i = 0 , 

2 sinh( / u) . 


x = 




sin(fc) , 


+00 


k = sinh 2 (/j) ^ 


e n tanh[/i(n — a;)] 


cosh[/i(n + 1 — x)] cosh[/r(n — 1 — a;)] 
k 


a = 2[cosh(/x) cos(fc) H— sinh(/i) sin(fc)] 

h 

+ sinh 2 (/r) ^ e„(n - x) tanh[/i(n - a:)] 


cosh[/i(n + 1 — x)] cosh [/i(n — 1 — x)] 


+oo 


sinh(/i) cosh(yu) 


cosh[/x(n + 1 — x)] cosh[/z(n — 1 — x)] 


(17c) 


(17d) 


Eqs. (fI7|) are integrated numerically with the initial conditions (J9J) . The algorithm used is the 
4th-order Runge-Kutta with At = 0.001, and the summations are truncated to a finite window 
around the center of mass of the soliton. The values of the parameters /i, k, x, and a are in¬ 
serted in (J2J) • We refer to this solution as the semianalytical solution. We compare this solution 
with the numerical integration of Eq. (jHJ) (referred as the numerical solution), and the results 
are presented in Fig. [0 In Fig. [TJa) we compare the center of mass coordinate x of ((21) and (fill) 
found in the semianalytical calculation with the numerical solution, and small acceleration is 
presented for this type of soliton [see also Fig. |T](f)] . It is found that on the average the velocity 
of the semianalytical solution is somewhat larger than the one found numerically. In Fig. QJb) 
the velocity of the center of mass and the second moment are presented. The plots for the 
velocity are zoomed in Figs. (U(c)-(e). We note that for t < 100 there is excellent agreement 
between the numerical and the semianalytical results. At the time t > 100 the second moment 
increases rapidly, therefore the approximation (fTTjl is not justified anymore, and large deviation 
between the two solutions are shown. In Fig. Q](f), the semianalytical velocity for longer time 
and its linear fit are presented. 


B. The regime 1 

As a representative example in this regime, we study a soliton with fi = 1, moving in one 
realization of the random potential with W = 0.1. We solve numerically Eq. ([HD with the initial 
condition Q. In this case the initial width of the soliton is 7 sites (N w = 7). The results are 
presented in Fig. [2j We find that the center of mass x moves monotonically to the right till 
a time t — T c when oscillations start [see Fig. (2](a)]. From Fig. |2](b) we see that the velocity 
decreases monotonically for t < T c and oscillates for t > T c . The period of these oscillations 
decreases with time. From Fig. (21(c) we see that the mass decreases with time and in the Erst 
stage this decrease is rapid, therefore the approximation (ITTj) fails. Finally the mass approaches 
a nonvanishing constant. The interesting phenomenon we End is the trapping of the soliton for 
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t > T c as a result of randomness. The particle aspect of this dynamics will be discussed in the 
next section. 


C. The regime /i<Cl 


In this regime, the soliton has a larger width, and it is easier to lose its mass through ra¬ 
diation. With a limit of zero randomness, Gamier 0 shows that the soliton propagates with 
its mass decreasing to zero, and its velocity to a nonvanishing constant. In fact, the radiation 
induces a remarkable deformation on the soliton profile after some time of propagation if the 
randomness is not weak enough. In Fig. [3] we present an example with /i = 0.5 [where the 
initial width of the soliton is 11 sites (N w = 11)] and W = 0.1. The initial condition is (J9]) . 
From Fig. [3](a) we see that the soliton spreads, and radiates its mass over 300 sites in the time 
t = 1000. From Figs. EJ^b) and (c) we conclude that the mass and velocity decrease. Such de¬ 
crease of velocity contains some short time intervals where the velocity oscillates approximately 
near a constant. The similar time intervals have been observed in the work of Franzosi et.al. 


15], as they studied the mobile discrete breathers propagating on very weak backgrounds in 


the framework of discrete NLSE. Their time intervals appear to be much longer, with weaker 
velocity oscillation, since their background perturbation is very weak (of the level ICG 4 ~ ICG 3 ). 
However, we have not observed the trapping behavior in this regime, especially before the soli¬ 
ton undergoes a large deformation. 


III. PARTICLE APPROACH FOR SOLITON TRAPPING IN DISORDERED AL 
LATTICE 


In this section we will study the question: Can an AL soliton in a weak random potential 
be considered as particle? 

We will focus on the soliton trapping in the second regime where fi ~ 1, and try to give a 
particle description of the trapping behavior. We start from the momentum (|5]), and assume 
it still to be the momentum for the model (JH]) when the random potential is weak. Taking a 
derivative on both sides of Eq. ([5]) with respect to t, and substituting Eq. (JS]) into the result, 
we can obtain 

rip 

~fa = 2 J2 Re (VV0n+l)( £ n - £n+l) , (18) 


where Re means the real part. For derivation see Appendix B. With the assumption that the 
soliton is particle-like, Eq. (ITS]) can be viewed as the variation rate of its momentum. On the 
other hand, dP/dt by ([7]) can be also written as 


dP 

dt 


dM s 

dt 


v + Mg-r^r 


d 2 x 
dt 2 


(19) 


Notice that, due to the mass radiation, the term dM s /dt in Eq. (TTUj) can not be neglected, 
especially before soliton trapping. Therefore, we can write the randomness-generated force in 
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(20) 


two ways, one is directly 


F\ 


M, 


d 2 x 
dt 2 




using and (TTUl) we derive 


F 2 


+oo 

2 Re(^nC+l)( £ n -£n+0 

n=—oo 


dM s 

dt 


( 21 ) 


The test of the particle-like picture is performed by comparing the forces F\ and F 2 presented 
in Fig. O Excellent agreement is found. These results strongly support the description of 
solitons as particles. In this picture with the force F 2 we associate work done on the soliton 
that decreases its kinetic energy with an effective potential 


U{t) = U 0 - F 2 vdt 


( 22 ) 


'to 


where Uq is an parameter which can be viewed as the initial energy to be determined as the 
constant that leads the mean of U, over the time interval of trapping in simulations, to be 
zero, i.e., Uq =< f t>T F 2 vdt >. With the same data, we plot both of this effective potential 
U(t) and soliton velocity in Fig. 0(e). It shows that the first reflection (T c ~ 1100), with the 
soliton velocity changing its sign, occurs at a peak position of the effective potential. Also other 
changes in the direction of motion of the soliton take place at maxima of U(t) as can be seen 
from Fig. 0(e). This is a direct result of (l22|i since 


dU(t) 

dt 


-F 2 v , 


therefore = 0 implies either F 2 


0 or v = 0. 


(23) 


IV. SCALING OF THE TRAPPING TIME T c 


In this section we demonstrate that there exists a scaling relation between the time T c when 
the trapping starts and the random potential strength W. In Figs. 0(a) and (b) T c is plotted 
as a function of W. It is found that 

T c ~ W~ v , (24) 

with rj = 2.32 ± 0.41. For each W G [0.06,0.1], we average the function U(t) of (j22| ) over 6 
different realizations to derive U aV e{t), and plot it till the minimum value of T c of 6 realizations 
in Fig. 0(c). From Fig. 0(d), we see that U ave is related to W and t via the combination tW 2 . 
This suggests the scaling relation 

U ave « T{tW 2 ) , (25) 

where T is the scaling function. If trapping starts at the same value of T, one finds 

T c oc W ~ 2 . (26) 


7 





Here we want to give some comments on T c . In forming (124|) of this section, we use one 
realization of random numbers uniformly distributed in [—1,1], but multiplied by the strength 
Wj 2, as the random potential. Since change of T c is obvious for small variation of W [see 
Figs. 0(a) and (b)], it can generally reflect the scaling relation. One may arrange multiple 
realizations of random potential to compute the mean value of T c , as well as its variance 
possibly, in the statistic sense. 

On the other hand, if the randomness is very weak (W < 0.04), Eq. (1261) seems to be no 
longer valid since the soliton can propagate without reflection for very long time as U ave ap¬ 
proaches constant (long-time simulation reveals such feature). In this regime, before trapping, 
the soliton with its mass as almost a constant, loses its velocity very slowly during the very 
long time propagation in the random potential. In particular, T c can reach to about 10 5 for 
W = 0.02, and much longer for weaker randomness. We consider this behavior of the soliton 
to be similar as the moving breather in the discrete NLS lattice fisl . 16] under very weak noise. 


IV. SUMMARY AND CONCLUSIONS 


Dynamics of solitons in random potentials was studied in the framework of the Ablowitz- 
Ladik model |2], In particular we explored the question when can a soliton be considered as 
a particle and what are the conditions for trapping of solitons in a random potential. The 
behavior was classified into three regimes specified by Eq. (fT6j) . In the regime /i 1 for short 
times the approximation (TT71) holds. In particular yU changes some time resulting in the change 
of the soliton width. This destroys the semianalytic solution resulting of (TTTT) as is clear from 
Fig. \H An improved approach will be subject of future studies. For /i <C 1 the soliton spreads 
very quickly and the potential picture is not appropriate. 

The most interesting regime is when fi ^ 1. The most interesting phenomenon is that the 
soliton is trapped and moves as a particle with varying mass. The equality of Fj = F 2 that is 
demonstrated in Fig. [4] is a strong evidence for the particle nature. The velocity changes its 
direction at some maxima of the potential (l22j) as can be seen from Fig. 0(e) and as expected 
from (l23jh Better understanding of the potential U{t) and its relation to the average of the 
random potential over the profile of the soliton will be left for future studies. Finally we found 
that T c , the first time when the velocity changes direction, scales with the strength of the ran¬ 
dom potential according to (T25li and the potential is scaling function of time and the strength 
of random potential (l24|h This may signal the existence of an underlaying statistical theory 
that should be explored in the future. 
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Appendix A 


Refs. [7[ show that for an AL model with a perturbation term 

.dfa 


dt 


— —fyn -1 + f/Vt+l) (1 + \lpn | ) + Rn , 


(A.l) 


the soliton parameters in (J2]) in the adiabatic approximation satisfy the following evolution 
equations 


+OO 


/i = sin. 


h(» 


cosh[p(n — x)]Im(r n ) 


^ cosh[/i(n + 1 — x)] cosh[/i(n — 1 — x)] 


2 sinh(u) sinh(u) 

x = -— sm(k) H-— 




+oo 

E 


+oo 


k = sinh(/i) 


(n — x) cosh[p(n — x)]Im(r n ) 
p ^ cosh [n(n + 1 — x)] cosh [fi(n — 1 — x)] 

sinh[//(n — x)]Re(r n ) 


^ cosh[/r(n + 1 — a;)] cosh[/x(n — 1 — x)] 
k 

a = 2[cosh(/i) cos(fc) 4— sinh(p) sin(fc)] 

v-' (n — x) sinh[/x(n — x)]Re(r n ) 

+ smh (/ij y _ 

n=— oo 


cosh[/i(n 4-1 — x)] cosh[/i (n — 1 — x)] 


+oo 


COS. 


k 


Mao 

—c 
+c 

E 


cosh[p(n — x)]Re(r„ 


_ cosh[//(n 4-1 — x)] cosh[/i(n — 1 — x)] 
sinh(/i) (n — x) cosh[/z(n — x)]Im(r n ) 


li cosh[/x(n + 1 — x)] cosh[/i(n — 1 — x)] 

where r n = R n exip[—ik(n — x) — ia]. For Eq. (JH]) with the solution form 

_ e n sinh(/i) 
n cosh[/i(n — x)] 

Substituting (A.3) into (A.2) we obtain Eqs. (1171) . 


we have 


(A.2a) 
(A.2b) 
(A.2c) 


(A.2d) 


(A.3) 


Appendix B 

In this appendix, we will show the derivation of Eq. (ITS]) . Take a derivative on both sides 
of Eq. (J5J) with respect to t, one can obtain 


dP 

dt 


= i 


+oo 

; E 


dfa ,* , , d^* n+l 


dt 


fan+1 + fa 


dt 


dfa, , * d ^ n+ 1 

'fa+1 - 


dt 


dt 


(B.l) 
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With Eq. OSJ), we derive the following sets 


d'lb 

-QjT = *[(^n-l + '0«+l)(l + V’nV’n) - eA] , 

= -*[«-! + C+lX 1 + iC) - £n<] , 

^' = X(Xn + V’n+2)(1 + V’n+iV’n+i) “ ^n+lXn+l] , 

'—jjf 1 = + C + 2 )( 1 + i+lC+l) - £ n+lC+l] • 

Substituting (B.2) into (B.l), after simplification, we can obtain 


dP 

dt 


+oo 

Y MVw/’n +1 + V’nV’n+l) ~ ^n+1 (V’nV’n+l + V’nV’n+l)] 

n=—oo 

+oo 

2 y Re (ic+i)( e n- e «+i) • 

n=—oo 


Thus Eq. (fTg|) is derived. 


(B.2a) 

(B.2b) 

(B.2c) 

(B.2d) 
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Figure captions and figures 
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Figure 1(b) 
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Figure 1(c) 



Figure 1(d) 
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Figure 1(e) 
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FIGURES 1: AL soliton acceleration by one realization of random potential [/j, = 3, v(t = 
0) = 1, and W = 0.1]. (a) Comparison of the center of mass between the semianalytical 

and numerical results, (b) Comparison of the semianalytical and numerical velocities, and the 
second moment m 2 ■ (c)-(e) Zoomed views of the velocity comparison in three different time 
intervals, (f) Semianalytical velocity by solving Eqs. (fT7|h and the dashed line is the linear fit 
of the data. 13 








































































































































































Figure 2(a) 



Figure 2(b) 



Figure 2(c) 

FIGURES 2: AL soliton trapping by one realization of random potential [// = 1, v(t = 0) = 1, 
and W = 0.1]. (a) Center of mass x [Eq. (fill] ] as function of time t. (b) The soliton velocity v 
[Eq. (USD] as function of time t. (c) The mass M s [Eq. ([TUP ] as function of time t. 
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Figure 3(a) 




FIGURES 3: AL soliton propagation in one realization of random potential [// = 0.5, v(t = 
0) = 1, and W = 0.1]. (a) The soliton profiles at t — 0 and t = 1000 (note that the site 
coordinate is fixed on the center of mass), (b) The mass M s [Eq. ([TO]) ] as function of time t. 
(c) The soliton velocity v [Eq. (1T3|) ] as function of time t. 
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FIGURES 4: (a) Comparison of the two forces F x [Eq. ODD ] and F 2 [Eq. Oil) ] with the 

soliton parameters and random potential the same as those in Figs. E](a) and (b). (b)-(d) 

Zoomed views of the panel (a) for three different time intervals, (e) The effective potential U 
and soliton velocity v after the soliton’s first reflection (T c ~ 1100). 
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Figure 5(a) 




Figure 5(b) 




Figure 5(c) Figure 5(d) 

FIGURES 5: (a) T c as function of W. The random potential are consisted of one realization 

of random numbers in [—1,1], multiplied by different strength W/ 2. (b) Linear fit of the data 
(logarithm forms of variables) in (a), (c) U ave as function of t and W. (d) U ave as a function of 
the scaling variable tW 2 . 
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